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PREFACE 

The  purpose  of  this  paper  is  to  show  that  the  turbulent  field  in 
a  variable-density  shear  layer  is  equivalent  to  one  of  a  conatant- 
density  layer  L  which  the  lateral  dimension  becomes  a  random  func¬ 
tion  of  time-  The  results  of  this  study  can  be  applied  to  the  inter 
pretation  of  laboratory  and  range  data  of  compressible  wakes-  This 
study  is  par'  of  RAND's  work  for  ARPA  on  the  basic  properties  of  re¬ 
entry  wakes*  The  author  is  Professor  of  Aerospace  Engineering  at  the 
University  of  Southern  California  and  a  Consultant  to  The  RAND  Corpo¬ 
ration- 
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ABSTRACT 


The  paper  is  concerned  with  the  calculation  of  the  mean  flow 
field  of  free  turbulent  layers  of  variable  density.  It  is  shown  that 
if  the  velocity  distribution  in  a  particular  constant-density  flow  is 
known,  it  is  possible  to  obtain  the  corresponding  variable-density  ve¬ 
locity  field  without  the  introduction  of  a  compressible  turbulent 
("eddy")  viscosity.  This  is  accomplished  by  a  Dorodnitsyn-Howar th 
type  of  transformation  applied  to  the  time-dependent  rather  than  to 
the  mean  equations  of  motion,  as  was  done  in  the  past.  When  the  trans¬ 
formed  equations  are  averaged,  using  Reynolds'  method,  the  incompres¬ 
sible  turbulent  equations  for  the  mean  flow  are  obtained.  These  equa¬ 
tions  can  then  be  handled  by  conventional  methods.  It  is  shown  that 
che  predictions  obtained  by  this  procedure  agree  well  with  experimen¬ 
tal  results. 
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LIST  OF  SYMBOLS 


C 

P 

D 

d 

‘(f) 

*(^) 

h 

L 


l 


M 

P 

<*k 

R 

T 

T_ 

T 

C 

U,V,W 

U|V,W 

X.Y.Z 


x,y,« 


specific  heat 

drag  or  thrust  force 

diameter  (of  a  jet  or  of  a  body) 

self -preserving  velocity  function  In  a  variable-density  field 

self-preservlnp,  cemperature  function 

enthalpy 

characteristic  length  of  the  velocity  field  In  a  constant- 
density  flow 

characteristic  length  of  the  velocity  field  In  a  variable- 
density  flow 

characteristic  lengch  of  the  temperature  field  in  a  variable- 
density  flow 
Mach  number 
pressure 

heat-flux  vector 

gas  constant 
temperature 
total  temperature 

time 

velocity  components  la  constant-density  flow 
velocity  components  In  variable-density  flow 
coordinates  In  the  constant-density  field 
coordinates  In  the  variable-density  field 
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F  ■  mass  concentration 
V  "  ratio  of  specific  heats 
9  ■  momentum  diameter 

0T  ■  energy  thickness 

p  *  density 


•  shearing  stress  tensor 

Subscripts 

c  -  along  axis 

o  -  characteristic  value  (u  *  u  -  u  ) 


1  •  jet  exit 

■  «  conditions  outside  of  shear  layer 


I.  Introduction 


The  mean  velocity  field  In  constant-density,  free  turbulent  shear 
layers  can  be  quite  satisfactorily  predicted  using  a  phenomenological 
approach  (especially  if  the  effect  of  Intermittency  is  included).  Such 
an  approach  assumes  a  turbulent  exchange  coefficient  (or  turbulent  vis¬ 
cosity)  that  is  constant  across  the  shear  layer-  In  a  variable-density 
layer,  however,  where  the  Reynolds  stresses  depend  on  the  density,  such 
an  assumption  la  clearly  unsatisfactory-  This  fact  constitutes  one  of 
the  main  difficulties  in  formulating  this  problem-  The  present  paper 
suggests  a  method  that  overcomes  this  difficulty- 

In  the  past,  an  appropriate  transformation  has  been  applied,  with 
12  3  4 

reasonable  success,  *  ’  ’  to  laminar,  compressible  shear  flows,  relat¬ 
ing  them  to  the  incompressible  case-  This  approach  has  also  been  used 
for  turbulent  flows.  In  fact,  there  is  a  large  body  of  literature  in 
which  compressible,  turbulent-flow  problems  are  treated  by  the  adoption 
of  a  Dorodnltayn  type  of  scaling  --as  used  In  the  laminar  case  -- 
and  of  some  rather  arbitrary  assumptions  concerning  the  compressible, 


turbulent- transport  terms*  The  first  of  these  attempts  was  made  by 
Mager.^  The  only  systematic  attempt  to  find  a  suitable  transformation 
was  made  by  Coies^  for  the  turbulent  boundary- layer  problem.  Inciden¬ 
tally,  in  almost  all  cases  examined  by  these  authors,  attention  was 
focused  on  an  incomplete  form  of  the  equations  of  the  mean  motion. 

In  formulating  the  problem  at  hand,  it  is  appropriate  to  inquire 
first  whether  in  fact  the  application  of  a  relatively  simple  kinematic 
transformation  of  the  Dorodnltsyn  type  is  reasonable  for  turbulent  flows. 
The  following  observations  are  pertinent  with  regard  to  this  point: 

(1)  Although  reliable  experiments  in  variable-density  flows  are 
still  few  in  number,  there  is  reasonable  evidence  that  the  structure 
of  the  turbulent  vortlclty  field  is  not  altered,  significantly  in  the 
presence  of  density  or  temperature  fluctuations.  This  suggests  that 
intersetion  between  the  vorticity  and  entropy  (or  concentration)  modes 
is  probably  not  strong,  even  in  flows  of  moderate  Mach  numbers.  Indeed, 
Chu  and  Kovaaznay7  have  shown  theoretically  that  i  a  homogeneous  field 
these  Interactions  are  of  second  order. 

(2)  The  mean  conservation  equations  indicate  that  the  momentum 
and  energy  equations  are  coupled  primarily  through  the  spatial  and 
temporal  variations  of  the  density  and  the  transport  properties  of  the 

g 

gas.  For  the  case  of  free  shear  layers,  in  which  the  direct  effects 
of  the  Reynolds  and  Prandtl  numbers  are  negligible,  at  least  for  moder¬ 
ate  Mach  numbers,  the  main  c.  ipllng  occurs  through  the  density  varia¬ 
tion  only. 

These  remarks  suggest  that  an  attempt  to  soak  a  kinematic  trans¬ 
formation  that  would  decouple  the  momentum  and  energy  equations  is  a 
reasonable  approach  to  the  problem. 
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The  present  paper  first  rederives  the  mem  momentum  and  energy 
equations  for  a  variable -density  turbulent  flow  in  a  form  most  conve¬ 
nient  for  comparison  to  the  constant -density  case*  In  this  derivation 
the  turbulent  transport  terms  are  generated  from  the  nonlinear  convec¬ 
tion  terms,  using  a  time-averaging  procedure.  These  nonlinear  terms, 
furthermore,  are  known  to  be  amenable  to  the  Dorodnitsyn-Howarth  type 
of  transformation  use.,  in  laminar  flows.  It  is  clear  that  past  diffi¬ 
culties  in  applying  such  a  transformation  to  turbulent  flows  arose  from 
the  time-averaged  value  of  these  terms-  This  suggests  an  alternative 
approach,  in  which  the  Dorodnitsyn-Howarth-Moore  transformation  is  ap¬ 
plied  to  the  time-dependent  equations  of  motion  and  the  time-averaging 
procedure  is  subsequently  carried  out  with  the  introduction  of  certain 
approximations,  a  method  that  allows  the  recovery  of  the  incompressible 
Reynolds  equations.  This  result  implies  that  if  the  mean  velocity  field 
in  a  constant-deuaity  flow  is  known,  the  corresponding  variable-density 
field  can  be  calculated  without  the  introduction  of  a  hypothesis  con¬ 
cerning  the  compressible  turbulent  viscosity.  We  Invoke  only  the  clas¬ 
sical  analogy  between  the  turbulent  momentum  and  heat  or  mass  transfer- 
Using  this  procedure,  calculations  made  for  several  shear  flows  indi¬ 
cate  good  agreement  with  experimental  observations. 

II.  Preliminary  Remarks 

The  conservation  equations  for  a  viscous,  compressible  fluid  may 

9 

be  written  in  the  following  form: 


» 


at  ax 


3pu^  dpu^u^ 
3t  + 


st  5*i 


3p  brtk 
dxt  *\ 


.  i£+u  *l_  +  t 

3c  j  3xj  Jk  3^  3x^ 


We  next  decompose  the  various  flow  quantities  into  a  mean  and  a  fluc¬ 
tuating  component,  following  Reynolds*  We  will  assume  that  the  pres¬ 
sure-fluctuation  level  p '/p  is  small  compared  to  that  of  the  velocity* 

Such  an  assumption  is  reasonable  in  flows  in  which  the  amplitude  of 
the  velocity  fluctuations  is  small  compared  to  the  velocity  of  sound; 
this  condition  is  satisfied  in  flows  of  moderate  Mach  number*  A  more 
detailed  diacusalon  of  this  point  is  given  by  Kistler**0 

In  introducing  the  velocity  and  temperature  perturbations  we  choose 
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a  maaa  weighted  average,  following  the  suggestion  of  Morkovin  and  Favre- 


Thus , 


u  -  u  +  u ' ,  p  -  p  +  p  ' 


where 


=  *  f 


and  P  is  the  conventional  time  average*  It  follows  that 


pu 7  “  0  and  P  “  0 
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Similarly,  we  write 


T  -  ?  +  t' 


(7) 


where 


T  = 


£T 

P 


(8) 


Substituting  Eqs .  (4)  and  (7)  into  the  perfect  gas  law 


p  *  RpT 


(9) 


(we  restrict  the  problem  to  this  case)  one  obtains,  after  averaging, 

P  *  RPT  (10) 


which  retains  the  form  of  the  perfect  gas  law.  Substituting  the  per- 
turbations  into  the  conservation  relations,  averaging,  and  using  bound¬ 
ary-layer  approximations,  one  obtains  the  following  equations: 


|£H+  pi  -  o 

ox  Oy 


(ID 


.  3u 

+  Pv  5— 
dy 


Spu  'v  ' 
3y 


(12) 


pu 


af  .  ~~  ax 

+  pv 

3x  3y 


(13) 


(The  molecular  viscosity  and  heat  conductivity  terms,  as  well  as  the 
mean  pressure  gradients,  have  been  neglected.) 

The  advantages  of  using  mass  weighted  averages  can  be  readfly 
seen.  First,  the  continuity  equation  does  not  have  a  source  term;  sec¬ 
ond,  the  Reynolds  and  energy  equations  retain  their  forms  corresponding 
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to  the  incompressible  esse-  It  la  further  seen  that  the  momentum  and 
energy  equations  are  coupled  only  through  the  density-  Thus,  a  coordi¬ 
nate  transformation  Involving  the  density  immediately  suggests  itself- 
In  fact,  since  the  convective  terms  have  the  same  form  as  those  for 
laminar  flows,  many  authors  have  applied  to  the  former  a  transformation 
used  for  the  laminar  case,  as  discussed  in  the  Introduction-  The  tur¬ 
bulent  transport  terms,  however  (since  they  contain  a  variable  density), 
cause  a  basic  difficulty  in  applying  the  transformation. 

The  main  purpose  of  this  paper  is  to  devise  a  method  to  overcome 
this  difficulty-  Our  method  takes  advantage  of  the  fact  that  the  tur¬ 
bulent  transport  qusntitiaa  arise  from  the  nonlinear  convective  terms, 
which  are  amenable  to  the  Dorodnitsyn  type  of  coordinate  transformation. 
This  fact  immediately  suggests  that  the  transformation  should  be  applied 
not  to  the  time-averaged  Reynolds  equations,  but  to  the  nonstationary 
conservation  Eqs-  (1)  and  (2)-  The  perturbations  are  subsequently  in¬ 
troduced  into  the  transformed  aquations  and  the  everaglng  procedure  is 
carried  out-  The  turbulant- transport  terms  generated  by  the  averaging 
process  are  found  to  be  Indeed  Independent  of  the  local  density- 

For  subsequent  calculations  it  is  useful  to  Introduce  some  well- 
known  integral  relations  of  the  above  equations.  By  integrating  Eq- 
(11)  with  respect  to  y,  one  obtains  the  momentum  integral  relation 

~  j**  pu(u  -  tO  dy  -  0  (14) 

where  U  is  the  velocity  outside  of  the  shear  layer.  Integrating  with 
respect  to  x,  we  obtain 
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pu(U  -  a)  dy  ■  const. 


(15) 


where  D  corresponds  to  the  wake  drag  on  a  body  for  a  wake  flow  or  to 
the  momentum  flux  through  an  orifice  for  a  jet  flow. 

A  similar  relation  is  obtained  from  the  thermal  energy  equation 

J  pu(T  -  T^)  dy  “  const.  ■  E  (16) 

where  C^E  is  the  excess  enthalpy  flux  in  the  wake  and  that  at  the  ori¬ 
fice  exit  for  a  jet  flow.  Since  a  constant  pressure  field  is  assumed, 
Eq.  (16)  can  be  rewritten  in  the  following  form  using  relation  (10): 

-  ».r 

oo  -m 

and  therefore 

P)  dy  -  E  (17) 


Another  useful  integral  relation  may  be  obtained  for  the  mean  ki¬ 
netic  energy  flux.  If  the  momentum  equation  (12)  is  multi.  '  <*d  by  u 
and  integrated  with  respect  to  y,  one  obtains 

J  pG(G  -  UJ2  dy  -  2  J  (U^  -  u)  dy  (18) 

•  00  •  OD 


III.  The  Formul  Transformation 

As  discussed  in  the  Introduction,  we  are  seeking  a  transformation 
that  would  eliminate  the  density  variation  in  the  conservation  equations 
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lor  mass  and  momentum.  To  this  end  we  Introduce  the  following  trans¬ 
formed  coordinates: 


X  -  x 


pjf(x,y,z,t)  -  |  p(x,y,z,t)  dy 


Z  -  z 


iSext ,  we  choose  velocity  components  in  the  transformed  plane  (say, 
the  incompressible  plane)  chat  satisfy  the  zero-divergence  conditions. 
This  may  be  accomplished  by  f'  e  following  relations 


P  V  -  PV  +  p  (~  +  W  —■  +  K 
*  \dt  ox  iz  / 


W  “  w 


Substituting  those  expressions  into  the  continuity  equation  (1)  and 
the  X-cvuponent  cf  the  mooentuus  equation,  one  obtain* 


M  a.  5V  .  3W  _  _ 


dU  ,  du2  .  3UV  ,  3W  . 

St  *  sr +  w *  *r  *  vl“tou#  •1CTrA3 


New  we  introduce  the  following  perturbations: 


U(X,Y,Z,t)  -  C(X.t)  +  u'fX,Y,Z,t] 
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V(X,Y,Z,t)  -  V(X,Y)  +  V  '[X, Y,X, t] 

(28) 

W(X,Y,Z,t)  -  W'[X,Y,Z,t] 

(29) 

Y(x,y,z,t)  -  Y(x,y)  +  Y'(x,y,z,t) 

(30) 

Here  the  bars  denote  time  averages, 


whence  by  definition 


U7»77«W7-77  =  0  (31) 

When  expressions  (27)  to  (30)  are  substituted  into  Eqs.  (25)  and  (26) 
and  averaging  is  carried  out,  terms  of  the  form  dt7dY  will  arise,  where 
$  denotes  a  velocity  component  or  products  thereof-  Since  both  $  and 
Y  are  random  functions  of  time,  the  mean  value  of  their  derivative 
shown  above  is  difficult  to  determine-  Lacking  adequate  statistical 
Information  concerning  these  functions,  we  shall  assume  that 


d$  d*  a? 

dY  **  dY  “  dY 


(32) 


To  clarify  the  Implication  of  such  an  approximation,  consider  the 
transformed  form  of  $,  say,  ^(XjV.Zjt);  then  for  the  two-dimensional 
case  using  Eq.  (20), 


— 'l  »  /dtp  dy\  m  dtp 
-dY /  Vdy  dY/  dy  p 


(33) 


In  a  single-component  system,  for  Instance,  p  /p  «  T/T  ;  hence 

CO  00 


&£  +  Stp'  T ' 

dy  T  dy  T 

40  00 


(34) 


4 


T 


This  procedure  was  suggested  to  the  author  by  Dr-  Y.  H.  Pao- 
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On  Che  ocher  hand, 


n  .  l. 

37  3y  37  3y 


1. 

oy  T 

OD 


(35) 


Thus  che  approximation  (32)  Implies  ChaC  Che  correlation  beCween  the 
fluctuations  of  che  temperature  and  velocity  gradients  la  neglected. 
Whether  chla  Is  justifiable  can  only  be  established  experimentally. 

Under  these  assumptions  the  time-averaged  forms  of  Eqs.  (25)  and 
(26)  become,  using  a  boundary- layer  type  of  approximation, 


3U  3V 
3X  i? 


3U2  SUV 

Sx~  &7 


o 


iPv7  +  viscous  terms 


(36) 

(37) 


This,  of  course,  is  the  incompressible  form  of  the  mean  flow  equations- 
It  is  seen,  therefore,  that  under  the  assumption  expressed  by  Eq-  (32), 
the  compressible  Reynolds  equation  may  be  transformed  into  the  Incom¬ 
pressible  form  for  free  shear  flows  with  a  constant  pressure  field- 
This  implies  that  the  transformation  (19)  to  (24)  provides  a  means  of 
obtaining  the  velocity  distribution  in  a  variable-density  shear  flow. 
According  to  Eq-  (22)  the  axial  component  of  che  velocity  profile  is 
in  fact  identical  to  that  of  the  corresponding  constant-density  flow, 
provided  that  the  lateral  coordinates  are  appropriately  stretched  in 
accordance  with  Eq-  (20).  To  carry  out  che  coordinate  transformation, 
on  the  other  hand,  we  need  to  know  the  density  (or  temperature)  dlstri- 
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button-  The  latter  may  be  obtained  by  assuming  the  classical  analogy 
between  momentum  and  mass  (or  heat)  transfer,  thereby  introducing  a 
turbulent  Prandtl  number  or  a  ratio  between  the  characteristic  length 
scales  of  the  two  processes*  The  only  undetermined  constant  in  the 
foregoing  method  is  connected  with  the  virtual  origin  of  the  shear  flow. 
Conditions  preceding  the  development  of  turbulent  shear  flows  are  gener- 
ally  complicated,  especially  if  a  laminar-to-turbulent  transition  is  in¬ 
volved.  However,  this  problem  is  beyond  the  scope  of  the  present  paper- 

Since  many  applications  Involve  an  axisymraetric  geometry,  the  fore¬ 
going  analysis  has  been  carried  out  for  this  case  also.  Although  the 
general  conclusions  are  the  same,  several  additional  points  should  be 
noted - 

The  coordinate  transformation  appropriate  to  an  axisymmetrlc  geom¬ 
etry  is 


X  -  x 


(19') 


I  P(x,y, i.  . t)y  dy 


(20) 


Y 


(21') 


where  y  and  t  are  the  radial  and  azimuthal  coordinates,  respectively. 
The  velocity  components  in  the  transformed  plane  will  have  the  follow 
lng  form; 


U  -  u 


(22  ') 


p  VY 
0» 


pvy  +  p  Y 

OD 


(£♦ 


dY 

dx 


+ 


W  dY\ 
Y  &*/ 


(23') 


y 


(24  ') 


Substituting  these  expressions  into  the  appropriate  forms  of  the  non- 
staclonary  compressible-continuity  and  axial -momentum  equations,  and 
following  the  same  steps  and  assumptions  described  above,  one  recovers 
the  Incompressible  form  of  the  conservation  equations. 

IV.  Applications 

In  Section  III  it  was  shown  that  by  use  of  the  transformations 
(19)  to  (24)  the  flow  field  in  a  variable-density  turbulent  shear  flow 
can  be  calculated,  provided  the  velocity  distribution  in  the  constant- 
density  case  is  known  and  an  assumption  is  made  concerning  the  turbu¬ 
lent  transport  of  heat  or  mass.  This  assumption  can  take  the  form 
either  of  an  empirical  value  of  the  turbulent  Prandtl  (or  Schmidt) 
number,  or  of  a  ratio  of  two  length  scales,  f  and  L^,  corresponding 
to  the  characteristic  scales  for  the  momentum  and  heat  (or  maaa)  trans¬ 
fer,  respectively-  These  conclusions  will  now  be  applied  to  two  types 
of  shear  flows,  the  axlsymmetrlc  jet  and  the  wake,  and  the  results 
compared  with  experimental  evidence. 

It  is  well  known  that  for  incompressible  shear  flows,  simple  inte¬ 
gral  techniques  can  predict  (except  for  an  undetermined  constant)  the 
spreading  rate  of  the  turbulent  zone  and  the  streamwlae  variation  of 
the  velocity  field  without  the  Introduction  of  any  phenomenological 
turbulence  theories.  This  approach  will  be  used  here. 

The  Circular  Jet 

Experiments  indicate  chat  both  Che  mean  velocity  and  mean  temper¬ 


ature  in  a  variable-density  jat  are  self-preaerving  except  near  the 
jet  exit  (see  Ref.  20,  p.  170).  Thus,  one  may  write 
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(38) 


1  -  *.♦*„«(£) 
T 


(39) 


Here  u  is  the  characteristic  (or  center)  velocity  and  T  s  T  -  T 
c  J  o  c  tj 

(the  difference  between  the  temperature  at  the  jet's  center  and  the 
reference  temperature  just  outside  of  the  turbulent  zone).  The  char¬ 
acteristic  length  l  is  defined  by  the  following  integral: 

h  *  f  ‘  Gf)  1  '  1  <40) 

o 


The  relation  between  £  and  will  be  discussed  later. 

The  momentum  integral  corresponding  to  (15)  for  a  circular  jet 
has  the  following  form: 

rm  ~  _2  2 

2n  J  pu  y  dy  *  P1°1A  ”  const-  (Al) 

0  O  ' 


The  subscript  1  corresponds  to  the  known  conditions  at  the  jet  exit, 
and  A  is  the  Jet  ares.  Applying  the  transfc  .ration  (20)  and  (22),  Eq. 
(Al)  becomes 


2np  U^L2 

ao  c 


*2  ®  l  4 


(A2) 


or 


U2 

il2 

U1 


^i  di  i_ 
8  L2 


_  e 


(43) 
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where 


e 


2 


(44) 


Note  that  the  momentum  diameter  9  is  the  proper  reference  length  of  the 

12 

problem,  as  pointed  ouc  by  Kleinstein-  In  Eq>  (42)  we  used  the  self- 
preserving  condition 

V  •  3  P  (*)  (45) 

where  L  is  the  characteristic  length  of  the  constant-density  jet  and, 
again, 

X1  K  f  F  ©t  4  “  1  (46) 

o 

Similarly, 

X2  .  f  P!  |  %  (47) 

0 

The  transformed  kinetic-energy  Integral  (18)  becomes,  for  a  cir¬ 
cular  Jet, 


dY 


(48) 


Assuming  now  that  the  Reynolds  stresses  also  follow  a  salf- preserving 


(49) 


distribution,  so  that 
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Equation  (AS)  simplifies  to 


;jx  uVi, 

dX  c  3 


-  2i?Uc 


(50) 


where 


S 


'3H 


*"<>  ic  ■  r«.)'4  (5i> 


With  Eq.  (A3),  Eq.  (50)  further  simplifies  to 


i-l!k 

-2  dX 
c 


'Jh  1 

■  I3  u? 


(52) 


which  becomes,  after  integration, 


U 


51  ‘  *  ^ 

C  ■» 


X  -  X 


(53) 


Substituting  Eq-  (53)  into  Eq-  (A3),  we  obtain 


L 

e 


2  XG  X  -  Xo 

2ri — 


(5A) 


From  measurements  of  constant-density  jets,  one  finds  that 


l&  t 

2  — -  -  0-052 

X3 


(55) 


4  13 

This  value  corresponds  to  one  given  by  Kinze,  based  on  the  exper- 

1A 

ieents  of  Van  der  Hegge  Zijnen,  if  the  fact  that  p./p  -  .91  for  these 
experiments  is  taken  into  account-  It  is  also  consistent  with  Corrsin 
and  Uberoi's  measurements,  where  P^/PB  “  -95. ^ 


-Iflm  mk  imibkilkmtlkeliUmkt;ijb(  asUtfJiUeUjk  ..^nUL;  Ui;j 
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The  value  of  Xq  depends  on  Che  flow  conditions  at  the  nozzle  exit,  as 
discussed  in  Section  III.  In  low-speed  flows,  however,  this  value  is 
small  compared  to  the  value  of  X  at  which  self-preservation  begins,  and 
can  generally  be  neglected.  A  good  approximation  to  the  values  of  the 
definite  integrals  and  may  be  obtained  by  assuming  a  Gaussian 
distribution  for  F(Y/L);  then  I^  *  1/2  and  I^  *  1/3.  Such  an  assump¬ 
tion  is  of  course  unnecessary  if  one  solves  the  equations  of  motion  di¬ 
rectly  using  a  phenomenological  approach.  Thus,  the  solution  for  the 
Incompressible  problem  has  the  form 

U.  X  -  X 

~  -  0.052  — g — 2.  (56) 

C 

X  -  X 

I  -  0.073  -y-S  (57) 

for  values  of  X  sufficiently  large  that  self-preservation  holds.  Using 
now  the  transformation  (19)  and  (22),  we  can  rewrite  Eq.  (56)  for  a  var¬ 
iable-density  Jet  as 


U,  x  -  x 

y  “  0-052  -g  (58) 

c 

The  characteristic  length  l  is  related  to  L  through  the  coordinate 
stretching  given  by  Eq*  (20).  This,  in  turn,  depends  on  the  density 
variation.  To  obtain  the  density  (or  temperature)  distribution,  we 
need  some  knowledge  of  the  turbulent  heat-transfer  mechanism.  If  we 
assume  Chat  the  mechanism  is  the  same  as  that  of  the  momentum  transfer, 
Chen  the  governing  Eqs.  (12)  and  (13)  have  the  same  form  and  g(y/l)  ■ 
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f(y //)  for  similar  boundary  conditions.  Experiments  with  conventional 

averaging  methods  show,  however,  that  the  characteristic  length  f.i 

the  normalized  temperature  distribution  is  larger  than  that  for  the 

2  2 

velocity.  Thus,  letting  a  B  L  /l^, 


m  ■ « © 


Assuming  a  Gaussian  distribution,  this  becomes 


From  the  energy  Integral  (16)  applied  to  the  axisymmetrlc  case  we  ob- 


00 

J  pu(T  -  T^ly  dy  -  plUl^Ti  *  Tro)A  “  const. 


or,  using  Eqs •  (38),  (39),  and  (60), 


u  T 


p  U.  T  .!  p  y  y  T  p  8  .2 

00  i  oo  o  m  0 


eT  T  -  T 
T  _  I  « 


where  0T  is  the  "energy  thickness"  of  the  Jet.  Applying  the  transfor¬ 
mation  (24)  to  the  integral  in  Eq.  (62),  we  obtain 


f  *  -  A.  f  ^  H  * 


Therefore  Eq.  (62)  becomes 


A  similar  expression  can  be  obtained  for  the  total  temperature  decay 
and  concentration  decay,  since  both  are  governed  by  an  Integral  equa¬ 
tion  similar  to  Eq.  (17).  Thus 


9.7  (O  +  1) 


(67) 


and  for  the  mass  concentration  at  the  jet's  center,  we  have 

rc  -  9.7  (o  +  i)  rt  j4'"x .  (68) 

o 

where  is  the  mass  concentration  at  the  orifice  exit  (the  concentra¬ 
tion  is  taken  as  zero  in  the  ambient  field). 

Finally,  in  order  to  obtain  the  spreading  rate  of  the  Jet,  we  use 
the  alternate  form  of  the  energy  Integral  [Eq.  (17)]  applied  to  the 
axlsymnetrlcal  cate: 
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or 


l/H  •  *i 


(70) 


But  the  value  of  Che  inCegrals  is  unicy  by  definition;  therefore 


L2  e2  L2  “c 


(71) 


Comparing  Eq.  (71)  Co  Eq.  (65),  we  obtain 


,2  i  * 

i_  1  «  * 

,  2  1  J+1T 

L  00 


(72) 


9.7 


*t  e 


e2  x  *  x0 


(73) 


or,  using  Bq.  (57), 


l 

e 


2 

2 


5*3  X  10° 


+  5.2  X  10"2 


®T  X  *  Xo 


(74) 


It  is  seen  in  Eq.  (74)  Chat  a  variable-density  jet  does  not  spread  lin¬ 
early.  However,  unless  Che  "energy  thickness"  is  large,  Che  nonlinear 
Cera  is  negligible  in  che  self-preserving  region  of  Che  flow. 

Equations  (58),  (66),  and  (72)  completely  determine  the  flow  fluid- 
The  only  "adjustable"  constant  is  xQ,  which  could  be  estimated  by  con¬ 
sidering  the  mixing  sons  and  Che  poCentlal  cone  near  Che  noszle  exit- 
Instead,  however,  we  shall  choose  Che  value  of  xo>  using  Eq.  (72),  chat 
gives  the  best  agreement  with  experimental  results.  Equations  (58)  and 
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(66)  can  Chen  be  unambiguously  compared  to  measurements*  In  all  cal¬ 
culations  we  use  a  m  0.64,  the  value  obtained  from  the  measurements  of 
Corrsln  and  Uberoi-^ 

In  Fig-  1  the  measurements  of  several  investigators  are  shown- 

The  supersonic  jet  experiments  of  Eggers^  and  Johannesan*7  extend 

well  beyond  the  self-preserving  region  of  the  jet,  and  Che  virtual 

18 

origin  was  easily  determined-  Only  two  sets  of  Warren's  data  are 
reproduced;  his  ocher  subsonic  measurements  show  similar  trends,  where¬ 
as  his  supersonic  ones  do  not  extend  far  enough  In  the  self-preserving 
region  to  be  usable-  The  selected  values  of  xq/6  are  Indicated  on  the 
figure;  the  origin  of  the  abscissa  was  shifted  for  clarity. 

Figure  2  shows  the  center  velocity  decay-  The  solid  line  corre¬ 
sponds  to  Eq.  (61).  It  is  seen  that  the  agreement  with  the  experi¬ 
ments  is  quite  satisfactory- 

In  Fig-  3  the  center- temperature  defect  is  plotted  against  the 

2  2 

axial-distance  parameter  8  x  -  x  / 9_9 -  The  solid  line  corresponds  to 

o  T 

Eqs -  (66),  (67),  and  (68);  in  all  cases  o  "  0-64-  Also  shown  are  the 

static  temperature  measurements  of  Corrsln  and  Uberol,  ^  the  total 

18  19 

temperature  data  of  Warren  and  the  results  of  Keagy  and  Weller 

obtained  in  a  helium  jet-  (Eggers  and  Johannesen  did  not  take  tem¬ 
perature  data.)  Although  these  experiments  are  comparatively  crude 
(the  ratio  of  probe  to  jet  width  is  large),  they  are  of  interest  be¬ 
cause  of  the  large  density  ratio  p  /p. -  The  agreement  between  the 

co  X 

analysis  and  the  measurements  is  again  satisfactory- 
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Fig.  1 — Jet  width 
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Fig.  2 — Center  velocity  In  the  circular  jet 


Fkj.3— Center  temperature  in  the  circular  jet 
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The  Axisymmetrlc  Wake 


It  1b  well  known  chat  for  wake- like  flows,  self-preservation  oc¬ 
curs  only  if  the  characteristic  velocity  uq  (the  centerline  velocity 
defect)  is  small  compared  to  the  free-stream  velocity.  Thus, 


u  »  u  -  u  f 
«  O 


for  —  «  1 


Under  this  assumption  the  temperature  distribution  is  also  self- 
similar,  having  the  form  given  by  Eq.  (39). 

20 

In  the  constant-density  case,  Townsend  has  already  shown  chat 

the  velocity  defect  and  wake  spreading  are  given  by  the  following 

—  2 

equations,  neglecting  terms  of  the  order  (U  / U  )  : 

O  oo 


v.  -  m  w 

»  -  fef 


Here  the  reference  length  6  is  defined  as 


9*  .  V 
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Cq  being  the  wake  drag  coefficient  of  the  body.  The  integrals  1^,  1^, 

and  I  ere  given  by  Eqs.  (40),  (47),  and  (51). 

g 

There  are  very  few  experimental  results  available  for  the  incom¬ 
pressible,  axlsymmetrlc  wake*  It  is  thus  difficult  to  obtain  a  reliable 
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value  for  the  constants  In  Eqs*  (76)  and  (77).  One  observation,  how- 

2 1 

ever,  should  be  made  using  the  results-  According  to  Reichardt, 
these  constants  differ  with  the  shape  of  the  body  considered.  For  in¬ 
stance,  the  spreading  rate  of  a  sphere  is  higher  than  that  of  a  cylin¬ 
drical  body  (with  its  axis  in  the  flow  direction).  One  should  there¬ 
fore  be  careful  in  comparing  constant-  and  variable-density  wake  flows 

For  the  case  of  a  cylindrical  body,  Hall  and  Hislop22  obtained 
20 

experimentally  the  value 


2 


14.1 


(79) 


Unfortunately,  since  their  measurements  have  been  carried  out  only  to 
about  seventeen  diameters  downstream  of  the  body,  where  self-preserva¬ 
tion  is  not  yet  completely  achieved,  the  above  value  might  be  somewhat 
high.  Using  the  value  tn  Eq.  (79)  and  with  ^  *  1,  Eqs .  (76)  and  (77) 
become 


(80) 


(81) 


For  the  variable-density  wake,  Eq.  (80)  retains  its  form,  and  tne 
decay  of  the  characteristic  temperature  can  again  be  calculated  from 
the  energy  integral*  Accordingly, 


2n  [  pu(T  -  TJy  dy 


E  e  2n0tp  U  X 

X  <*  T  OD 


(82) 
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By  using  arguments  similar  Co  chose  for  Che  jet  problem,  we  may  write 


¥  -  T 
00 

'm  — 

T 
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(83) 


so  fchaC  Eq-  (82)  becomes 


f  L2 
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(U 


0  F )F°  y  d y 
o  L  L 


e£u  x 

X  a>  oo 


(84) 


This  simplifies  co 


82  2 

_T  0^ 

e2  l2 
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(85) 


1  - 


a  +  1  u 


where  again  a  Gaussian  distribution  has  been  assumed  for  the  velocity 
profile- 

Xo  express  the  "energy  diameter"  in  terms  of  free-stream  condi¬ 
tions,  one  assumes  that  all  Che  kinetic  energy  loss  in  Che  wake,  as 
Indicated  by  the  momentum  defect,  is  transformed  into  heat-  Thus, 
using  Eqs-  (15),  (78),  and  (82),  we  obcaln 

DU  -  E  (86) 

00 


or 


2np  IfV  -  2ttP  U  C  X  62 

00  oo  03  OD  p  ID  X 


(87) 


Xhls  gives 
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i  -  (V  -  DM2 


Substituting  this  expression  into  Eq.  (85)  and  using  Eqs-  (80) 
snd  (81),  one  obtains 


-  2.8o(v  -  1)M‘ 


I  /X  -  x  v ' 


x  -  x.\-2/3 


,  2.8  a  /X  '  o\ 

1  *  rr^(-r”) 


Prom  the  axlsynmetric  form  of  Eq.  (17)  we  again  obtain  the  width 


of  the  jet: 


2nT.u„  f  (p.  -  dy  -  2ttVooso  J  fy  dy 


+  2rrT 


■=  r 

•  o  J 


fpy  dy  -  2tt€,p  U  T 

T  •  ">  oo 


But  since 


P.  -  P  ’  P  x2  ** 


Equation  (90)  becomes 


5s4l+!..2  . 

at  L  \T  1  +0^L 

®  <E  ■ 


or,  using  Eq.  (85),  one  recovers  Eq.  (72)  for  the  Jet;  that  is, 


1  To 
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(93) 


Iil^i  te-JMiitr  li-cliiims. 


With  Eqs •  (80),  (81),  and  (85)  this  finally  becomes 
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These  results  are  compared  to  the  measurements  of  Demetrlades  In 
Fig-  4.  Since  Demetrlades  presents  his  data  In  terms  of  x  -  x^,  Eqs- 
(80),  (89),  and  (93)  can  be  directly  related  to  the  measurements- 
Again  a  -  0-64  was  used  for  the  computations  -  The  agreement  between 
the  present  analysis  and  the  experiments  is  acceptable- 

V.  Conclusion 

The  analysis  carried  out  In  this  paper  Indicates  that  with  the 
use  of  a  suitable  transformation  it  Is  possible  to  predict  with  rea¬ 
sonable  accuracy  the  mean  velocity  and  density  (or  temperature)  field 
In  a  variable-density  shear  layer,  once  the  velocity  distribution  in 
the  corresponding  constant  density  layer  is  known- 

Two  explicit  results  follow  from  the  transformation: 

(1)  The  characteristic  velocity  of  a  self-similar  shear  layer 
is  a  function  of  the  nondlmensional  coordinate  x/B  only,  Irrespective 
of  the  density  variation  In  the  flow-  Here  6  Is  a  reference  length 
defined  by  the  momentum  integral  equation- 

(2)  The  ratio  of  the  width  of  a  variable-density  layer  to  that 
of  the  corresponding  constant-density  layer  can  be  expressed  explicit¬ 
ly  in  terms  of  the  local  maximum  temperature  (or  concentration)  differ¬ 
ence-  For  axially  symmetric  geometries,  this  has  the  form 
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1  + 


a  +  1  T 

CO 


There  are  several  obvious  llmicacions  of  the  method; 

(1)  A  constant  pressure  field  is  assumed.  This  assumption  might 
be  relaxed  --at  least  for  certain  types  of  pressure  gradients  --if 

a  more  sophisticated  transformation  were  used. 

(2)  The  early  stages  of  the  shear- layer  development  cannot  be 
treated,  especially  if  a  laminar- to- turbulent  transition  occurs.  One 
is  then  forced  to  assume  a  "virtual  origin"  for  the  shear  layer. 

On  the  other  hand,  the  proposed  method  opens  new  possibilities 
for  the  analysis  of  compressible,  turbulent  boundary  layers. 
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